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ABSTRACT 


A  now  technique  for  the  correction  of  atmospheric  distortion  in  a 
telescope  is  presented.  Most  of  this  distortion  arises  from  a  random 
phase  variation  in  the  incoming  light  across  the  te’escope  aperture. 
This  variation  limits  the  resolving  power  of  even  large  telescopes  to 
about  one  second  of  arc.  If  one  defines  the  sharpness  value  of  the 
images  in  a  suitable  way,  this  sharpness  is  maximized  only  when  the 
phase  distortion  of  the  incoming  light  is  zero.  The  report  presents 
computer  simulations  of  a  simple  feedback  system  in  which  flexible 
optical  elements,  adjusted  to  maximize  the  sharpness,  correct  most  of 
the  atmospheric  distortion.  Photon  statistics  set  the  limiting  magni¬ 
tude  of  the  object  for  which  a  practical  feedback  system  can  work.  One 
should  bo  able  to  resolve  details  in  a  7th-magnitude  object  to  better 
than  0.1  second  of  arc.  The  system  can  be  conveniently  employed  within 
existing  telescopes. 
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The  work  reported  here  was  jointly  supported  by  the  Defense 
Advanced  Research  Projects  Agency,  the  National  Aeronautics  and  Space 
Administration,  and  the  Atomic  Energy  Commission.  The  results  will 
also  appear  in  a  report  issued  by  the  Lawrence  Berkeley  Laboratories 
(LBL-2492)  ,  and  will  appear  in  the  Journal  of  the  Optical  Society 
of  America. 
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Satisfactory  Definitions  of  Image -Plane  Sharpness  . 


I  INTRODUCTION 


The  resolving  power  of  optical  astronomical  telescopes  with 

aperture  greater  than  about  10  cm  has  been  limited  by  atmospheric 

distortion  rather  than  by  diffraction,  ever  since  the  first  astro- 

1* 

nomical  telescopes  were  built  hundreds  of  years  ago.  Except  on 
unusually  calm  nights  at  the  best  observatory  locations,  atmospheric 
"seeing"  limits  the  resolving  angle  of  large  telescopes  to  about  1 
second  of  arc,  the  diffraction  limit  of  a  10-cm-diameter  telescope. 

This  limit  has  rendered  impossible  the  imaging  of  the  disks  of  nearby 
stars,  close  binary  systems,  and  the  detailed  structure  of  the  cores 
of  galaxies.  It  has  also  prevented  the  high-quality  imaging  ol  the 
planets  that  would  otherwise  be  achievable  with  existing  large  tele¬ 
scopes. 

Attempts  to  lessen  atmospheric  distortion  by  locating  the 
observatory  at  high  altitudes  have  been  only  partially  successful. 
Airborne  observatories  have  seeing  limitations  similar  to  those  on 
mountaintops,  attributed  to  boundary-layer  turbulence  close  to  the 

2  3 

airplane.  Balloon-borne  systems  avoid  atmospheric  distortion,  but 
the  difficulties  associated  with  balloon-flight  operations  make  this 
technique  impractical  for  an  extended  program  ot  astronomical  measure¬ 
ments.  Until  we  enter  a  new  era  of  large  telescopes  in  space,  atmos¬ 
pheric  distortion  will  remain  a  major  concern  of  astronomers. 

Michel  son  and  Pease  were  the  first  to  overcome  the  1  second-of-arc 
limit.4  Although  their  stellar  interferometer  was  incapable  of  complete 


*  References  are  listed  at  the  end  of  the  report. 
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image  restoration,  it  permitted  the  first  measurement';  of  stellar 
diameters,  close  binary  star  systems,  and  the  size  of  Jupiter's  moons. 

Considerable  improvements  on  interferometer  techniques  have  been  made 

5  6 

by  Hanbury  Brown  and  Twiss  and  by  Labeyrie  and  co-workers.  These 

interferometer  techniques  are  limited  by  the  fact  that  in  effect  they 
measure  only  the  intensities  of  the  spatial  Fourier  components  of  the 
object's  brightness  distribution.  Because  the  relative  phase  informa¬ 
tion  of  the  wave  components  is  lost,  only  the  auto-correlation  function 
of  the  image  can  be  constructed.  Unless  the  object  being  studied  has 

certain  symmetries  the  image  information  cannot  be  unfolded.  Multiple- 

7,8 

optical-interferometer  systems  overcome  this  limitation,  but, 

because  of  their  complexity,  they  have  not  yet  been  applied  to  astro- 

g 

nomical  objects.  Code  has  presented  a  review  of  recent  progress  with 
these  techniques.  In  principle,  if  one  knew  the  modulation  transfer 
function  of  the  atmosphere,  one  could  correct  the  distorted  image 
after  recording  it.  Considerable  progress  has  been  made  with  this 

g 

post-detection  approach,  but  because  the  modulation  transfer  function 
varies  with  position  and  with  time,  it  is  difficult  to  improve  resolu¬ 
tion  by  more  than  a  factor  of  2  for  astronomical  objects. 

An  alternative  to  post-detection  processing  is  real-time  compensa¬ 
tion  of  the  telescope  optics  system  in  order  to  cancel  out  the  phase 
distortion  introduced  by  the  atmosphere.  Several  workers  have  success¬ 
fully  used  devices  to  stabilize  the  motion  of  the  centroid  of  the 
image. 10  Unfortunately,  this  technique  is  incapable  of  providing  the 
substantial  improvement  necessary  to  yield  diffraction-limited  perform¬ 
ance  for  all  but  small  telescopes.  Babcock11  has  suggested  that  an 
active  corrector  plate  could  be  introduced  into  the  telescope  optics  to 
compensate  for  the  changing  atmospheric  phase  shifts.  In  order  to 
determine  tue  amount  of  correction  to  apply  to  a  given  region  of  the 
telescope  objective,  Babcock  suggested  performing  a  knife-edge  test  on 


2 


each  portion  of  the  objective,  using  an  unresolved  nearby  bright  star 
as  the  light  source.  Unfortunately,  this  technique  can  work  only  if  the 
light  from  the  nearby  star  has  experienced  essentially  the  same  atmos¬ 
pheric  phase-shift  perturbation  on  its  path  to  a  given  place  on  the 
telescope  objective  as  the  light  from  the  object  under  study.  Measure¬ 
ments  of  the  correlations  of  the  distortions  of  resolved  double 
12,13 


stars 


as  well  as  the  usual  theoretical  models  of  atmospheric  dis- 


14,15 

turbances  indicate  that  this  equal  perturbation  begins  to  fail  for 

objects  being  viewed  by  ground-based  telescopes  that  subtend  more  than 

a  few  seconds  of  arc.  The  angular  region  for  which  the  requirement 

holds  is  called  the  "isoplanatic  patch."  Since  the  atmospheric  distor- 

16 

tion  changes  with  a  time  constant  of  about  0.02  sec,  application  of 
Babcock's  scheme  to  astronomical  measurements  requires  that  a  bright 
unresolved  star  (typically  7th  magnitude  <•  v  brighter)  lie  within  a  few 
seconds  of  arc  of  the  object  being  studied.  Such  nearby  stars  are  rare. 
A  system  similar  to  that  proposed  by  Babcock  has  been  designed  by 
Perkin-Elmer  Corporation  for  a  different  problem:  in  situ  figuring  of 
a  light-weight  space  telescope  mirror.*7 


The  technique  we  introduce  in  this  report  was  conceived  in 

collaboration  with  F.  J.  Dyson  of  the  Institute  for  Advanced  Study, 

18 

Princeton.  Dyson  s  analysis  will  be  published  separately. 


Our  optical  system  requires  an  optical  element  that  can  shift  the 
effective  phase  of  the  light  across  the  aperture  of  the  telescope. 

This  could  be  accomplished  either  by  dividing  the  telescope  objective 
into  separately  movable  segments,  or,  more  conveniently,  by  introducing 
a  movable  optical  element  elsewhere  in  the  system.  Our  technique 
requires  no  unresolved  bright  star  nearby,  for  we  use  the  light  from 
the  object  under  observation  in  order  to  determine  the  required  phase 
shifts.  In  principle  our  technique  works  for  an  arbitrarily  complex 
object,  although  the  object  must  lie  within  the  isoplanatic  patch. 


3 


We  now  introduce  the  concept  of  image  "sharpness."  We  can  define 
this  sharpness  in  such  a  way  that  under  certain  conditions  the  value  of 
the  sharpness  for  an  atmospherically  degraded  image  is  always  less  than 
that  of  the  true  image.  One  such  sharpness  definition  is 

Si  =  /“dx dy  l“(x,y>,  (1) 

where  x,y  denote  coordinates  in  the  image  plane,  and  l(x,y)  is  the 
image  intensity.  Figure  1  is  a  schematic  illustration  of  our  real-time 
correction  scheme.  We  place  a  device  capable  of  measurement  of  the 
image  sharpness  in  the  image  plane  of  a  telescope.  The  disturbing 
atmosphere  degrades  the  image  sharpness  by  randomly  perturbing  the  phase 
and  amplitude  of  the  incoming  light.  Later  in  this  report  we  show  that 
the  amplitude  perturbations  are  unimportant  for  image  sharpness.  A 
phase  corrector  placed  in  the  optical  path  is  capable  of  restoring  the 
proper  phase.  Individual  segments  of  the  phase  corrector  are  adjusted 
in  turn  to  maximize  the  sharpness  of  the  image,  as  measured  by  the 
device  placed  in  the  image  plane.  If  the  entire  telescope  objective 
can  be  corrected  during  the  period  of  time  in  which  the  atmospheric 
distortion  holds  constant,  the  resulting  image  will  be  nearly  diffrac¬ 
tion-limited. 

In  the  remainder  of  this  report  we  discuss  several  alternative 
definitions  of  sharpness.  In  the  Appendix  we  present  proofs  that  some 
of  these  sharpness  functions  reach  their  maximum  only  for  a  properly 
restored  image,  one  from  which  the  atmospherically  introduced  phase 
variation  has  been  removed.  We  discuss  various  means  of  choosing  the 
most  practical  sharpness  functions.  We  then  present  computer  simula¬ 
tions  of  realistic  telescope  systems,  showing  that  our  feedback  scheme 
properly  restores  images  that  have  been  degraded  by  simulated  atmos¬ 
pheric  phase  distortions.  Finally,  we  present  the  results  of 
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Introducing  statistical  fluctuations  into  the  sharpness-measuring 
device,  showing  that  such  fluctuations  set  the  ultimate  limitation  on 
object  brightness  for  which  our  technique  produces  diffraction-limited 
performance . 
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I 


I 

II  IMAGE- PLANE  SHARPNESS  FUNCTIONS 

A  good  definition  of  the  sharpness  S  of  an  image  is  one  for  which 
S  reaches  a  maximum  value  only  for  a  true  (undistorted)  image.  For  the 
pirposes  of  this  report  we  consider  only  "aperture-plane  distortion," 
defined  below.  Let  W(u,v)  be  the  instantaneous  wavefront  at  the  tele¬ 
scope  aperture  in  the  absence  of  any  distortion.  If  the  distorted 
wavefront  W’(u.v)  is  given  by 

W'(u.v)  =  W ( u , v )  •  J(u,v)  (2) 

then  we  have  aperture-plane  distortion.  Here  A(u,v)  is  any  complex 
function  of  the  position  (u,v)  in  the  plane  of  the  aperture.  Spherical 
aberration  can  be  expressed  in  this  way  and  is  therefore  an  example  of 
aperture  plane  distortion.  Under  certain  conditions,  atmospheric  dis¬ 
tortion  can  also  be  expressed  by  this  formula.  If  the  object  being 
viewed  by  a  ground-based  telescope  is  of  sufficiently  small  angular 
extent  that  light  from  all  parts  of  the  object  incident  on  the  same 
point  (u,v)  has  passed  through  essentially  the  same  perturbing  atmos¬ 
phere,  then  we  have  aperture-plane  distortion.  Such  an  object  is  said 
to  lie  within  a  single  "isoplanatic  patch." 

The  flexible  optical  element  of  Figure  1  allows  us  to  add  extra 
phase-shifts  to  the  Mu,v),  in  hopes  of  reducing  or  eliminating  the 
distortion.  We  want  to  define  the  sharpness  S  in  such  a  way  that  any 
Mu,v)  (which  now  includes  the  effects  of  correction  as  well  as  the 

original  distortion)  other  than  a  simple  translation  of  the  image 

[  ik(a+bu+cv)  "i 

Mu,v)  =  e  ,  where  a,  b,  and  c  are  constantsl  will  reduce 

the  value  of  S. 
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Table  1 


SATISFACTORY  DEFINITIONS  OF  IMAGE-PLANE  SHARPNESS 


Definition 


Computer 

Simulation 


Comments ^ 


S^  =  /dx  dy  I 

S  =  I  (x  ,  y  ) 

2  o  o 


Satisfactory 


Satisfactory 


Proven  (see  Appendix) . 

Makes  poor  use  of  photon¬ 
counting  statistics,  so 
satisfactory  only  for  bright 
objects . 


S^  =  / dx  dy  MI 


m+n  2 

S„  =  /dX  dy  !iiOi 

4  Sx%n 


Satisfactory  for 
M  =  round  hole, 

I  =  single  or 
multiple  stars 

Untried 


Proven  if  M  = 
tor ted  image. 


=  undis- 
Reduces  to 


S  if  M  =  *(x  ,y  ). 
2  o  o 


Proven . 


S  =/dx  dy  1 
5 

n>2 


Satisfactory  for 
n  =  2,  3,  j 


Proven  only  for  unresolved 


S  -  -  / dx  dy  I  r 

b 

2  2  2 
r  =  x  +y 


S  =  -/dx  dy  In  (I.’) 
7 


Satisfactory 


dx  dy  I  ln(l) 


Moment  of  inertia  function. 
Proven.  Reduces  ails  of 
light  distribution  but  gives 
poor  central  maximum.  A 
special  case  of  S^. 

Minimizes  the  "entropy"  of 
the  image.  Quantization  of 
I  (photons)  gives  meaning  to 


S  =  - /dx  dy  1 1 - 1  I 
8  1  o' 


Untried 


Defect  function. 
Proven . 


1 


In  preparing  this  report  for  publication  we  found  that  at  least  S| ,  ,  and 

S6  have  appeared  before  in  the  optics  literature,  although  we  had  not  seen 
them  when  we  did  this  work.  See,  for  example,  E.  H.  Linfoot. 

*Many  of  the  proofs  are  due,  in  whole  or  in  part,  to  Freeman  Dyson. 


There  are  many  definitions  of  image  sharpness  S  which  satisfy  this 
criterion.  An  example  (mentioned  in  the  introduction)  is 
S  =  /dx  dy  I^(x,y),  where  l(x,y)  is  the  intensity  at  a  point  (x,y)  in 
the  image  plane  of  the  telescope  system.  In  the  Appendix  we  prove  that 
for  monochromati c  light,  reaches  its  maximum  value  only  when  Avu.v) 
is  reduced  to  no  more  than  the  image  translation  described  above.  The 
proof  is  done  by  using  a  version  of  the  Fresnel-Kirchof f  equation  to 
calculate  the  intensity.  is  maximized  for  zero  distortion  irrespec¬ 

tive  of  the  object  intensity  distribution.  We  can  therefore  use  to 
provide  feedback  to  the  active  optical  elements  and  remove  the  distor¬ 
tion,  even  for  complex  objects  such  as  planets  and  galactic  cores  (if 
they  lie  within  an  isoplanatic  patch).  There  is  no  way  that  aperture 
distortion  can  increase  the  value  of  beyond  its  undistorted  value 
even  for  an  arbitrarily  complex  object. 

The  second  law  of  thermodynamics  suggests  another  definition  of 

sharpness.  Let  S,  =  I (x  ,y  ),  the  brightness  of  the  image  at  an 
2  o  o 

arbitrary  point  (x  ,y  )  on  the  image  plane.  For  objects  that  Lave  a 
o  o 

brightest  spot,  we  believe  S  will  be  maximized  when  the  image  of  that 

bright  spot  is  shifted  over  to  the  point  (x  ,y  )  and  also  the  distortion 

o  o 

is  eliminated.  We  have  not  obtained  an  analytic  proof. 

Yet  a  third  definition  of  sharpness  is  conceptually  similar  to  the 
above  two  definitions.  Let  S  =  /dx  dy  M(x,y) *  I (x,y) ,  where  M(x,y)  is 
the  transmittance  of  a  mask  placed  over  the  image.  In  the  case  in  which 
M  is  an  accurate  replica  of  the  true  undistorted  image,  S  is  a  good 

O 

definition  of  sharpness  that  reduces  to  above  when  the  distortion  is 

eliminated.  In  the  case  in  which  M  is  a  delta  function  at  (x  ,y  ),  S, 

o  o  e 

reduce=  to  S  above.  We  are  investigating  how  closely  M(x,y)  must  match 
2 

the  restored  image  in  order  for  S  to  be  a  good  definition  of  sharpness, 
and  what  sort  of  final  images  may  result  when  the  match  is  poor. 
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We  have  investigated  several  other  definitions  of  S.  In  Table  1  we 
list  tne  ones  we  believe  to  be  good  definitions  of  image  sharpness  in 
the  sense  that  they  reach  a  maximum  whf>n  image  distortion  (excluding  a 
simple  translation)  has  been  removed.  We  also  indicate  which  functions 
have  been  tested  in  our  computer  simulations  (to  be  described  in  Section 
IV),  and  which  functions  have  a  formal  proof  in  the  Appendix. 
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Ill  REAL-TIME  CORRECTION  OF  TELESCOPE  OPTICS 


In  the  previous  section  we  introduced  the  concept  of  image  sharp¬ 
ness  S,  where  S  reaches  a  maximum  value  only  when  the  aperture-plane 
distortion  has  been  eliminated.  In  this  section  we  describe  a  feedback 
scheme  in  which  small  optical— path— length  perturbations  are  introduced 
in  the  optical  path  as  in  Figure  1.  The  effect  of  these  perturbations 
on  S  is  detected,  and  the  appropriate  compensation  to  maximize  S  is 
determined  and  then  permanently  introduced  in  the  optical  path. 

For  the  case  of  astronomical  seeing,  we  can  make  several 

simplifications.  The  distortion  function  A(u,v)  typically  has  a 

coherence  length  of  10  cm,  corresponding  to  the  1  arc  sec  mentioned 

previously.  We  plan  to  use  discrete  compensating  segments,  each  one 

of  which  affects  approximately  10  cm  of  the  incoming  wave.  Therefore 

2 

a  telescope  of  aperture  area  A  cm  must  have  approximately  A/100 
compensating  elements.  Of  course  the  incident  wave's  distortion 
function  is  continuous  in  (u,v)  so  the  use  of  discrete  compensating 
segments  leads  to  a  residual  error  in  the  correction.  In  the  next 
section,  however,  we  show  that  this  residual  error  is  small. 

Another  simplification  results  because  "random  apodization"  of  an 

aperture  has  little  distorting  effect  on  an  image.  This  allows  us  to 
ikft(u,v) 

write  A(u,v)  =  e  where  6(u,v)  is  purely  real.  In  the  next 

section,  by  using  computer  simulation  we  show  that  ignoring  amplitude 
changes  in  the  incident  wave  due  to  the  disturbing  atmosphere  causes 
little  image  distortion.  It  is  not  difficult  to  see  why  this  is  so: 
the  central  maximum  of  the  image  of  a  point  object  has  all  of  the 
amplitudes  adding  up  in  phase,  whereas  at  the  first  minimum  these 


11 


amplitudes  all  cancel.  A  random  variation  of  amplitudes  at  the  tele¬ 
scope  aperture  somewhat  diminishes  the  intensity  at  the  central  maximum, 
corresponding  roughly  to  the  fraction  of  the  aperture  that  has  been 
blacked  out.  The  first  minimum  also  gains  some  intensity,  but  for  a 
telescope  whose  aperture  contains  many  10-cm  "cells"  the  gain  in 
intensity  is  small  compared  with  the  remaining  central  maximum.  On  the 
other  hand,  a  random  phase  variation  across  the  aperture  completely 
destroys  the  central  maximum.  In  tact,  the  average  intensity  expected 
at  the  original  site  of  the  first  minimum  and  the  original  site  of  the 
central  maximum  are  now  approximately  the  same. 

We  plan,  then,  to  use  a  compensating  system  that  alters  only  the 

phase  of  the  incoming  wave  and  not  the  amplitude.  The  adjustable 

11 

element  could  be  an  eidophor  system  as  proposed  by  Babcock,  a  set  of 

movable  mirror  segments  driven  by  piezoelectric  or  magnetostrictive 

transducers,  or  a  plate  of  variable  optical  thickness  (such  as  PLZT 
2() 

ceramic  materials  ). 

We  must  correct  the  telescope  objective  within  a  time  less  than 

the  characteristic  time  t  in  which  the  atmospheric  distortions  change. 

Measurements  of  the  frequency  spectrum  of  stellar  scintillations  and  of 

16 

stellar  speckle  patterns  indicate  that,  at  a  good  telescope  site,  we 

should  have  about  t  -  0.02  sec  to  make  our  corrections.  j  for  a  high- 

altitude  aircraft  should  be  about  ten  times  shorter.  Statistical 

fluctuations  in  the  number  of  photons  arriving  within  the  time  period 

set  the  minimum  photon  flux  of  an  object  whose  image  can  be  corrected 

by  this  technique.  Unless  we  have  at  least  one  detected  photon  for 

each  aperture  segment  in  the  time  t,  there  is  no  possibility  of  making 

2 

a  correction.  (For  t  =  0.02  sec,  and  100  cm  aperture  segments,  this 
argument  implies  that  the  object  under  observation  must  be  15th 
magnitude  or  brighter).  For  a  real  system,  considerably  more  than  one 
photon  per  segment  per  t  is  required.  We  will  now  make  an  order-of- 
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magnitude  estimate  for  the  minimum  flux  of  the  object  under  observation 
required  for  a  real  correction  system  to  work. 

Consider  that  n  settings  of  an  individual  correction  segment  are 
necessary  to  determine  the  proper  setting  of  that  segment  to  maximize 
S.  If  the  area  of  the  telescope  is  A,  and  the  effective  area  of  an 
individual  correction  segment  is  a,  then  the  number  of  individual  seg¬ 
ments  is  N  =  A/a.  If  n  settings  are  required  to  position  an  individual 
segment  (n  =  2,  typically),  the  time  available  for  each  setting  is 

t  =  T/nN  =  Ta/nA  .  (3) 

Let  the  flux  of  the  object  under  observation  be  B,  measured  in  terms  of 
2 

photons/cm  sec)  at  the  aperture.  Then  the  number  p  of  detected  photo¬ 
electrons  in  time  t  is 

p  =  TlBAt  =  TjBaT/n  (4) 

where  T|  is  the  photon  detection  efficiency.  Statistical  fluctuations 
dp  in  p  are  typically  ~  /p.  Our  computer  simulations  show  that  the 
fractional  change  in  a  sharpness  value  caused  by  varying  one  segment  is 
approximately  1/N,  although  a  poor  choice  of  the  sharpness  function  can 
cause  it  to  be  less  than  this.  In  order  to  make  a  statistically  signif¬ 
icant  measurement  of  the  effect  of  the  variation  we  must  have 

1/N  ^  dp/p  =  1/  /TlBai/n  .  (5) 

Substituting  N  =  A/a  we  find  that 

B  ^  nA  /a  Tjr  .  (6) 

We  can  combine  Eqs.  (4)  and  (6)  to  get,  simply, 

p  £  N2  .  (7) 
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In  the  next  section  we  verify  this  formula  empirically  with  computer 

4  2  2  2 

simulations.  For  A  =  10  cm  ,  a  =  10  cm  ,  T  =  0.25,  n  =  2,  and  t  =  0.02 

4  2 

sec,  we  find  B  M  x  10  photons/cm  sec.  Since  the  relationship 
between  astronomical  magnitude  m  and  B  (for  a  0.3  Urn  bandpass)  is* 


B 


(4  x 


6 

10  ) 


x  10 


-m/2.5 


(8) 


this  corresponds  to  objects  brighter  than  about  5th  magnitude.  This 
is,  of  course,  only  a  rough  estimate.  By  arranging  the  telescope 
geometry  in  the  form  of  an  annulus,  one  could  achieve  an  angular  reso¬ 
lution  of  about  0.02  arc  sec  for  this  example  (see  Section  V.). 


21 

♦Can  be  calculated  from  the  quantities  given  by  C.  W.  Allen. 
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IV  COMPUTER  SIMULATIONS 


We  have  simulated  telescope  systems  with  active  feedback,  using 
the  Fresnel-Kirchof f  equation  in  a  numerical  calculation  of  the  image 
intensity  distribution  I(x,y)  for  monochromatic  light.  In  the  cases 
where  we  formed  images  of  complex  objects,  we  assumed  that  the  object 
flux  distribution  I ( p, &)  was  spatially  incoherent.  To  simulate  the 
atmospheric  distortion  we  introduced  a  random-walk  phase  variation  into 
the  incident  wavefront  at  the  telescope  aperture.  The  characteristic 
distance  in  which  this  phase  variation  had  become  randomized  was  set  to 
be  10  cm.  Essentially  identical  results  were  given  by  an  alternative 
simulation  of  the  distortion,  in  which  the  incoming  wave  was  divided 
into  10- cm  portions  each  of  whose  phase  was  randomly  varied  between 
±n  radians. 

To  provide  the  simulation  of  actual  telescope  feedback  systems  we 
considered  the  aperture  of  the  telescope  to  be  broken  up  into  N 
separately  movable  segments.  We  introduced  a  small  phase  shift  in  the 
portion  of  the  incident  wavefront  corresponding  to  a  single  segment  and 
calculated  the  changed  value  of  the  sharpness  function  resulting  from 
the  perturbation.  The  measurement  of  sharpness  change  was  used  to 
drive  the  feedback  system,  determining  and  setting  a  permanent  phase 
correction  for  that  particular  aperture  segment. 

We  used  two  feedback  schemes.  In  the  first  we  introduced  one  and 
then  a  second  phase  addition  to  the  particular  segment,  in  order  to  get 
sharpness  values  corresponding  to  perturbations  of  +n/2  and  then  -n/2 
radians.  Letting  S(0)  be  the  original  sharpness  value,  and  S(+-)  ,  and 
S(-)  be  the  calculated  sharpness  values  for  the  shifted  segment  phase 


positions,  we  then  determined  the  best  setting  for  that  segment,  phase 
shift  cp,  using 

tan  cp  =  [SCO  -  S(-)]/lS(  +  )  f  S(-)  -  2S(0)  ]  .  (9) 

Equation  (9)  follows  from  the  fact  that  all  the  sharpness  functions  we 
investigated  vary  nearly  sinusoidally  as  the  phase  shift  of  any  parti¬ 
cular  segment  varies  over  its  range  of  2 jt  radians. 

We  also  used  a  feedback  scheme  that  followed  the  gradient  of  the 

sharpness  function.  This  scheme  has  the  advantage  of  being  extremely 

simple  to  build  into  a  practical  telescope  system.  We  inti oduced  a 

small  perturbation  of  radians  into  an  individual  segment.  If  the 

sharpness  value  increased  we  left  the  perturbation  in,  and  mov^d  on  to 

the  next  segment.  If  the  sharpness  value  decreased,  we  removed  the 

perturbation,  placing  the  segment  at  -cp  instead,  before  moving  onto 

o 

the  next  segment.  Although  this  feedback  scheme  did  not  provide  the 

accurate  segment  setting  of  Eq.  (9)  ,  it  did  perform  nearly  as  well 

(for  cp  =  it/ 4)  when  allowed  to  cycle  through  the  N  mirrors  twice, 
o 

However,  most  of  the  calculations  we  shall  report  in  this  section  used 
the  technique  of  Eq.  (9). 

We  chose  N  =  25  aperture  segments  in  order  to  simulate  a  rela¬ 
tively  small  telescope  and  yet  be  able  to  achieve  a  resolution 
significantly  better  than  1  arc  sec.  We  considered  two  telescope 
geometries.  The  first  was  a  strip  or  "flatland"  mirror  that  we  took 
to  be  250  cm  long  by  10  cm  wide.  The  images  formed  by  such  a  mirror 
have  much  better  resolution  along  the  direction  of  the  strip,  so  we 
suppressed  the  other  dimension,  presenting  our  calculations  in  just 
a  one-dimensional  image  "plane."  The  second  geometry  was  a  ring 
mirror  with  an  outer  diameter  of  75  cm  and  an  inner  diameter  of  55  cm. 
For  this  telescope,  we  calculated  and  displayed  two-dimensional  images. 
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We  first  used  our  computer  program  to  confirm  our  claim  that  phase 
variations  in  the  atmospheric  disturbance  are  much  more  disruptive  to 
the  image  than  are  amplitude  variations.  Figure  2  shows  the  images  with 
the  strip  mirror  of  a  single  unresolved  star  for  (1)  no  atmospheric 
disturbance,  (2)  random  amplitude  variations  in  the  incident  wavelront, 
and  (3)  random  phase  variations  in  the  incident  wavefront.  It  is  clear 
from  Figure  2  that  the  amplitude  variations  merely  reduce  the  height  of 
the  central  maximum,  while  the  phase  variations  completely  disrupt  the 
central  maximum,  giving  a  speckle  pattern  similar  to  those  reported  in 
real  observations  by  Labeyrie.6  In  this  pattern,  the  size  of  the  over¬ 
all  image  corresponds  to  the  seeing  limit  (1  arc  sec) ,  whereas  uhe 
average  size  of  individual  speckles  is  close  to  the  diffraction  limit 
of  the  telescope. 

An  "iteration"  of  segment  settings  consists  of  a  single  pass 
through  the  N  segments,  setting  each  individual  phase  according  to 
Eq.  (9).  Figure  3  shows  how  the  atmospherically  perturbed  image  of  an 
unresolved  star  in  the  stri-  telescope  improved  as  each  of  the  25  seg¬ 
ments  was  adjusted  in  turn.  The  sharpness  function  chosen  for  this 
particular  calculation  was  from  Table  1.  An  image  < hat  is  essentially 
diffraction-limited  resulted  after  only  a  single  iteration.  In  the 
restored  image  approximately  half  of  the  original  flux  is  still  spread 
around  in  residual  speckles  whose  brightness  is  typically  1 0%  that  of 
the  newly  formed  central  maximum.  The  central  maximum  formed  at  the 
location  of  the  brightest  original  speckle  when  sharpness  functions  , 

S  or  S  were  used,  but,  as  expected,  it  formed  at  the  origin  or 
5  7 

(x  , y  )  when  we  used  function  S  ,  S  ,  or  S  .  Systems  employing  the 
o  o  2  3b 

latter  functions  use  the  phase-shifter  array  to  compensate  telescope 
tracking  error  as  well  as  the  atmospheric  perturbations.  Of  course, 
telescope  tracking  can  also  be  easily  accomplished  with  a  separate  feed¬ 
back  system.  Figure  4  shows  the  original  speckle  pattern  and  restored 
image  for  the  ring  mirror  geometry  and  sharpness  function  5,^  using  for 
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SECTION  A-A  INTENSITY 


(a)  NO  DISTURBING  ATMOSPHERE 


FIGURE  4  TWO-DIMENSIONAL  IMAGES  FROM  A  COMPUTER  SIMULATION 
OF  A  RING-SHAPED  TELESCOPE  APERTURE.  Sections  of  the 
contour  plots  are  shown  in  the  accompanying  histograms.  The 
contour  plots  (but  not  the  histograms)  have  been  smoothed  to 
average  our  unphysical  features  due  to  the  discrete  nature  of  the 
computer  calculations. 
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SECTION  B-B  INTENSITY 
—  percent 


SPECKLE  PATTERN  FOR  A  RANDOM  PHASE  DISTORTION  AT  THE 
APERTURE  PLANE 


FIGURE  4  (Continued) 
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SECTION  C-C  INTENSITY 


FIGURE  4  (Concluded) 


the  mask  M  a  hole  the  size  of  the  restored  diffraction  central  maximum. 
The  improvement  in  the  image  after  only  a  single  iteration  is  dramatic. 
Finally,  Figure  5  shows  the  original  speckle  pattern  and  restored 
image  for  three  stars  of  equal  flux  separated  by  angles  of  0.04  aid 
0.20  second  of  arc,  viewed  by  the  strip  mirror. 

The  residual  speckles  visible  in  the  restored  images  of  Figures  3 
through  5  were  not  reduced  further  by  subsequent  iterations  of  the 
mirror  system.  Such  a  residue  of  speckles  in  the  restored  image  is  not 
unexpected,  both  in  the  simulations  and  in  practice,  since  the  perturbed 
wave  varies  continuously  over  the  aperture,  whereas  we  have  applied  the 
corrections  to  N  discrete  segments  in  the  movable  optical  element. 

Figure  6  shows  the  initial  phase  distribution  of  the  atmospherically 
disturbed  light  from  an  unresolved  star  at  the  telescope  aperture, 
tr  ether  with  the  N  values  of  cp  calculated  for  this  case  according  to 
Eq.  (9).  It  is  clear  from  Figure  6  that  our  applied  correction  is  only 
approximate.  Therefore  perfect  restoration  of  the  proper  straight  wave- 
front  is  impossible  with  only  25  discrete  elements,  and  our  applied 
correction  cannot  completely  restore  the  image.  We  allowed  this  mis¬ 
match  in  our  simulations  because  we  felt  that  such  a  discrete  compensa¬ 
tor  might  be  easiest  to  build,  using  separately  movable  mirrors  or 
transparent  phase  shifters.  A  single  continuous  mirror,  deformed  by 
discrete  transducers,  may  have  smaller  residual  speckles  than  indicated 
in  Figures  3  through  5,  since  it  exploits  the  continuity  that  must  be 
present  in  any  real  incoming  wavefront. 

Most  of  the  sharpness  functions  we  tried  from  Table  1  converged  on 

the  same  solution  in  one  (or  at  most  two)  iterations.  An  exception, 

however,  was  S  =  f  r^I,  which  weighted  the  far-away  speckles  so  heavily 
6 

that  a  central  maximum  could  not  form.  Figure  7  shows  the  result  of  two 

iterations  with  S  .  Reduction  of  light  in  the  "tails"  of  an  image  may 
6 


23 


SPECKLE  PATTERN 


PHASE  ANGLE, 


SEGMENT  PHASE 
SETTINGS  AFTER 
ONE  ITERATION 


10  15  20 

APERTURE  CORRECTION  SEGMENT  NUMBER 


FIGURE  6  PHASE  ANGLE  vs,  POSITION  FOR  THE  RING-TELESCOPE 
GEOMETRY.  Position  is  indicated  by  numbering  the  25 
correction  segments  around  the  ring.  The  original  disturbed 
phase  ov  the  incident  wave  (smooth  curve)  and  the  segment 
setting  resuming  for  a  single  iteration  (histogram)  are  shown. 
The  image  resulting  from  this  wavefront  was  shown  in 
Figure  4.  Segments  1  and  2  were  off  by  27r  and  have  been 
redrawn  (dashed)  to  show  another  equivalent  position. 


SPECKLE  PATTERN 


be  valuable  for  some  applications,  but  in  most  cases  the  central  maximum 

formed  by  S  is  too  poor  to  be  useful. 

6 


The  easiest  sharpness  function  to  implement  on  a  real  telescope  is 

S  ,  since  a  device  to  make  this  measurement  consists  merely  of  a  single 

photomultiplier  with  a  mask  of  the  desired  transmittance  M(x,y)  placed 

in  front  of  it  at  the  image  plane.  Our  proof  in  the  Appendix  applies 

only  in  the  case  in  which  M(x,y)  corresponds  to  the  true  restored  image 

I  (x,y),  but  we  have  found  in  the  computer  simulations  that  S  still 
o  3 

works  very  well  when  M  is  only  approximately  equal  to  I  .  In  particular, 


we  tried  an  M  that  was  a  round  hole  in  an  otherwise  opaque  mask:  the 
diameter  of  the  hole  was  chosen  to  correspond  to  the  size  of  the  diffrac¬ 
tion-limited  image.  This  definition  of  S  gave  results  as  satisfactory 

o 

as  those  we  achieved  with  S  or  S  .  All  sharpness  functions  we  tried 

1  5 

gave  as  good  an  image  restoration  for  complex  objects  (double  stars  and 
triple  stars)  as  for  a  single  unresolved  star.  In  practice,  the  choice 
of  sharpness  function  depends  on  both  ease  of  implementation  and  the 
type  of  object  being  observed.  For  objects  such  as  binary  star  systems 
where  the  telescope  could  center  on  the  brighter  star  of  the  pair,  one 
would  probably  use  the  simple  Sr  system.  For  a  rapidly  tumbling  object 

w 

such  as  an  irregular  asteroid  the  sharpness  function  might  be  more 
appropriate.  For  imaging  the  disk  of  a  planet  a  function  related  to  S 

*3 

would  probably  be  best,  in  which  one  maximizes  the  light  through  a 
narrow  annular  slot  in  the  image  plane. 


We  used  our  computer  simulations  to  evaluate  the  effect  of  statis¬ 
tical  fluctuations  in  the  number  of  detected  photons  (photoelectrons). 

For  each  measurement  of  S  we  perturbed  the  number  of  photoelectrons  in 
each  image  bin  in  a  random  way,  using  a  Poisson  distribution.  Figure  8 
shows  the  resulting  image  of  the  triple  star  (each  star  of  9th  magnitude) 
after  one  iteration.  For  this  image  we  used  the  sharpness  function  , 
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a  strip  mirror  of  25  segments  (each  100  cm  ) ,  a  quantum  efficiency 
J\  =  0.25,  and  an  atmospheric  coherence  time  t  =  0.02  sec.  The  average 
number  of  photoelectrons  available  for  each  calculation  of  S  was  1000. 
When  the  number  of  photoelectrons  was  reduced  by  a  factor  of  three,  the 
image  was  severely  degraded. 

We  define  the  image  quality  Q  to  be  S/S^,  where  S  is  the  sharpness 

value  achieved  in  the  presence  of  statistical  fluctuations,  and  S  is 

0 

the  value  achieved  if  the  statistical  fluctuations  were  absent.  A  plot 
of  Q  versus  number  of  photoelectrons  is  shown  in  Figure  9  for  our  simu¬ 
lations  of  an  annular  telescope  consisting  of  25  segments.  Again, 
roughly  1000  photoelectrons  per  segment  setting  are  needed  in  order  to 

make  a  significant  improvement  in  the  image  sharpness.  Equation  (7) 

2  2 

would  have  predicted  that  N  =  25  =  625  photoelectrons  are  required. 
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fluctuations.  The  sharpness  function  used  was  S3  =  f  Ml,  where  M  =  a  round  hole 
the  same  size  as  the  central  diffraction  pattern  of  a  point  star.  The  number  of 
photoelectrons  was  defined  to  be  0.25  times  the  number  of  photons  incident  on  the 
entire  telescope  in  one  setting  time  t.  Very  similar  results  were  found  for  Si  and 
for  the  strip  telescope  geometiy. 


V  DISCUSSION  AND  CONCLUSIONS 


We  have  shown  in  the  preceding  sections  that  a  basically  simple 
system  of  active  telescope  optics,  using  a  sharpness-measuring  device 
at  the  telescope  image  plane  to  feed  back  to  a  flexible  optical  element, 
should  eliminate  most  of  the  image  degradation  caused  by  atmospheric 
distortion.  A  great  variety  of  image-sharpness  definitions  can  in 
principle  drive  the  feedback  system.  The  bandwidth  of  the  flexible 
optical  element,  set  by  the  characteristic  time  t  in  which  the  atmos¬ 
pheric  disturbance  changes,  appears  to  be  within  present-day  technical 
reach.  The  angular  size  of  an  object  that  can  be  corrected  this  way  is 
limited  by  the  size  of  the  isoplanatic  patch.  The  angular  extent  of 
this  patch  depends  not  only  on  the  physical  size  of  the  disturbed  region, 

but  also  on  its  nearness  to  the  telescope.  For  ground-based  telescopes, 

12,13 

the  size  of  the  patch  is  only  a  few  seconds  of  arc,  whereas  for 

airborne  telescopes,  where  the  turbulent  layer  is  very  close  to  the 
telescope,  the  size  of  the  patch  should  extend  out  to  many  minutes  of 
arc.  This  would  be  suitable  for  imaging  large  objects  such  as  Jupiter. 

The  time  t  also  determines  the  minimum  brightness  of  an  object  that 
can  be  corrected  by  this  technique.  Equation  (6)  relates  the  object 
flux  B  and  t  to  the  particular  parameters  of  the  telescope  system: 


B  > 


2 

n  A 
3 

a  tT 


where 


2 

B  =  Object  flux  in  photons/cm  sec  at  telescope 

n  =  Number  of  segment  settings  needed  to  determine 
optimum  position  for  that  segment 
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A  =  Area  (cm  )  of  telescope  aperture 
2 

a  =s  Area  (cm  )  of  each  adjustable  segment 

T|  =  Quantum  efficiency  of  detector 

T  =  Coherence  time  of  the  atmospheric  distortion. 

2 

For  ground-based  telescopes,  a  is  typically  100  cm  and  r  ^0.02  sec. 
Using  Eq.  (8)  to  relate  B  to  stellar  magnitude  and  choosing  an  annular 
telescope  geometry  to  give  the  best  resolution  for  a  given  A  (the 
Rayleigh  criterion  for  an  annular  ring  gives  P  =  0.76X/D),  we  can  derive 
the  following  relationship  between  object  magnitude  and  the  best  resolu¬ 
tion  obtainable  assuming  telescope  parameters  matched  to  the  object: 


P  = 


-3  m/5 

10  \  10 


/  nr/n 


where 

fl  =  Best  resolution  (radians),  Rayleigh  criterion 
X  =  Wavelength  of  light,  cm 
m  =  Magnitude  of  object  under  observation. 

-5 

As  an  example,  take  m  =  7,  \  =  5  x  10  cm,  T  =  0.25,  t  =  0.02  sec,  a  = 

2  -7 

100  cm  ,  and  n  =  2  settings.  These  yield  a  R  of  3  x  10  rad  «0.06  arc 

sec.  This  resolution  would  be  achieved  with  an  annular  telescope  with 

a  diameter  of  %  130  cm  %  50  inches.  Of  course,  one  cannot  usually  match 

telescope  parameters  to  particular  objects.  We  present  the  calculation 

above  only  as  a  rough  estimate  of  the  resolution  obtainable  for  a  giver 

magnitude  with  our  system.  A  superior  system,  with  a  dimmer  limiting 

magnitude  for  a  given  resolution,  could  be  built  by  employing  multiple- 

* 

image  systems,  each  with  its  own  sharpness-measuring  feedback  scheme. 


Multiple-image  systems  have  been  suggested  to  us  by  Jack  Franck  and 
22 

Freeman  Dyson. 
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Note  that  the  m  that  appears  in  our  equations  is  the  magnitude  of  the 
entire  object  (assuming  it  is  all  within  an  isoplanatic  pa^ch),  not  the 
magnitude  of  the  feature  under  observation.  Thus,  for  example,  the 
light  trom  a  bright  star  can  be  used  for  image  sharpening,  even  if  the 
object  under  study  is  a  dim  companion. 

A  number  of  difficulties  may  arise  if  polychromatic  light  is  used. 
First,  the  sharpness  theorems  were  proven  in  the  appendix  just  for  mono¬ 
chromatic  light.  If  the  object  is  multicolored  it  is  possible  to 

r  2  rj 
increase  -J  I  with  a  suitable  aperture-plane  distortion.-  However, 

L  23 

it  can  be  shown  that  the  theorem  for  is  still  correct  for  poly¬ 
chromatic  light  as  long  as  the  object  is  uniformly  colored — i.e.,  each 
point  on  the  object  emits  light  with  the  same  spectrum.  Second,  the 
lateral  color  aberration  pr<xluced  by  the  atmosphere  will  disperse  the 
image  unless  one  is  viewing  close  to  the  zenith.  Hill  and  Zanon!^  have 
proposed  a  telescope  modification  that  adequately  corrects  for  this. 
Third,  there  is  the  possibility  that  an  aperture  segment  in  our  system 
reaches  a  setting  different  from  the  true  setting  by  2n  radians.  Figure 
6  shows  two  segments  for  which  this  happened.  For  monochromatic  light 
such  a  shift  makes  no  difference,  but  for  polychromatic  light  it  causes 
some  image  degradation.  Of  course  a  narrow-bandpass  filter  would  solve 
these  problems,  but  at  a  cost  in  least  correctable  magnitude.  For  the 
third  problem,  we  expect  that  the  telescope  feedback  system  will  "lock 
in  on  the  correct  phase-distortion  distributions  as  they  drift  across 
the  aperture  of  the  telescope.  In  any  case,  if  one  wanted,  the  tele¬ 
scope  could  have  a  "search  mode"  in  which  2n  phase  shifts  are  also  tried 
on  the  individual  movable  segments  to  see  if  any  improvement  lies  in 
that  direction. 

The  image-sharpening  technique  discussed  here  can  be  used  for  more 
than  just  the  elimination  of  atmospheric  distortion.  Since  maximizing 
of  S  automatically  refigures  in  the  telescope  optics,  sharpening  could 


be  used  to  refigure  lightweight  optics  either  on  the  ground  or  in  space. 
Image  sharpening  can  also  be  applied  to  radio-interferometric  measurements 
where,  for  long-baseline  multielement  systems,  the  relative  phases  of 
the  receivers  are  unknown.  These  phases  can  be  determined  in  the  compu¬ 
ter  analysis  by  calculating  the  values  that  maximize  the  sharpness  of 
the  radio  image. 

Our  research  group  at  Berkeley  is  currently  designing  a  compensat¬ 
ing  system  that  could  be  attached  near  the  focal  plane  of  existing 
optical  telescopes. 
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FOR  THE  RESTORED  IMAGE 
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PROOFS  THAT  SHARPNESS  FUNCTIONS  REACH  MAXIMA 
FOR  THE  RESTORED  IMAGE 


Consider  an  image  plane  (x,y)  where  we  observe  an  intensity  distri¬ 
bution  I(x,y)  of  focused  light  from  a  distant  object  of  small  angular 
extent.  The  focusing  is  provided  by  a  telescope  objective  placed  a 
distance  f  away.  Let  (u,v)  be  the  coordinates  at  the  aperture  (objective) 
plane,  and  consider  the  objective  and  the  image  planes  both  to  be  per¬ 
pendicular  to  tiie  line  connecting  the  center  of  the  objective  and  the 
center  of  the  light  distribution  of  the  distant  object.  The  intensity 
I ( x , y )  is  then  given  by  a  modified  E  resnel-Ki rchoff  integral  over  the 
surface  of  the  telescope  objective: 


I  (x,y)  =  y”d0  d0  0(0  ,0) 


/du  dv 


ik[  6  (u,v)  -t-(uxt-vy) 
e 


1  >(ud  +v0  )  I 


2 


( A- 1 ) 
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where  we  have  neglected  the  angular  obliquity  factors.  Here  k  is  the 
wavenumber  of  the  light  (which  we  take  to  be  monochromatic),  0(0,0)  is 
the  light  flux  from  angles  (0,0)  falling  on  the  telescope  objective, 
and  we  have  assumed  that  the  object  under  observation  is  radiating 
incoherently.  The  real  function  6(u,v)  includes  both  the  effects  of 
the  disturbing  atmosphere  and  of  the  correcting  optical  elements.  As 
we  have  stated  in  the  main  text,  we  have  ignored  the  fact  that  the 
atmosphere  modulates  the  amplitude  as  well  as  the  phase  of  the  incoming 
wave.  The  assumption  that  the  entire  object  falls  within  an  isoplanatic 
patch  is  equivalent  to  the  assumption  that  the  atmospheric  distortion 
can  be  written  in  the  form  Mu.v)  instead  of  the  more  general  form 
6  (u,  v,e  ,0) . 
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1.  Proof  that  Sharpness  Functions  and  S4  Reach  Absolute 
Maxima  for  a  Restored  Image 

Sj  is  a  special  case  of  the  more  funeral  sharpness  iunction  S^, 
defined  in  Table  1  of  the  main  text: 


S  =  f dx  dy 
4 


j!+m  , 

d  I  ( x ,  y ) 


m 


dx  dy 


(A-2) 


For  convenience,  we  rewrite  Eq.  (A-l)  eliminating  the  square  by  intio- 
ducing  the  dummy  variables  u'  and  v': 

I(x,y)  =  /d9  d0  0(6,0)  /du  du'  dv  dv ' 

ik[6(u,v)-6(u' ,v')+(u-u')x  f+(v-v')y  f] 

•  e 

^ik[ (u-u/)8+(v-v,)0 ]  (A_ 


We  now  substitute  Eq.  (A-3)  into  Eq.  (A-2),  perform  the  differentiation, 
and  write  out  the  square  by  introducing  new  dummy  variables: 


m 


S4  =  /dxdy  /d6i  d02  d0t  O(01,01)  0(e2.^2>  duj  duJ  du2  dU2  dVl  dVl  dV2  dV2 

2(£+m)  ,  l  .  ,  m  ,  l 

•  (k/f)  (u^u')  (vi>v1)  (u2'u2)  <v2_ 

ik[ (u^-u ' ) (x/f+0 1)+(vi-v')(y/f+0l) ] 


ik[  (u2-up  (x/f+e2)+(v2-v^)  (y/f+e2)  ] 


ik{5(ui,vi)-6(u',v')+6(u2>v2)-^(u^,vp} 


( A-4 ) 


The  integral  over  x  and  y  can  now  be  performed,  yielding  Dirac  delta- 

6D(urul+u2+u2>  and  6D(vrv;+vv2)-  The  integral  over 


funt  v  ions : 


these  delta-functions  can  be  performed  immediately,  eliminating  one  each 
of  the  u  and  v  variables.  Defining  new  variables  z  =  u^-u'  and 

w  =  v  -v ' ,  we  obtain 

1  1 


S  =  (?  )  /d0  d0  d0  d/»  0(8  ,0)0(6  ,0  )  /  du  dz  du  dv  dw  dvr 

4  IZIZIIZZ  X  z  1  ^ 

ik[(ft  -0  )z+(0  -0  )w] 

„  ,x  2(X+m-l)  2jC  2m  1  1  2  1^2 

•  (k  f)  z  w  e 


ik{  6  ( u  ,  v  )  -6  (u  -z ,  v  -w)  +*  ( u  ,  v  ) -6  ( u  +z ,  v  +w)  j 

XX  XX  m  fa  £* 


(A-5) 


Since  the  term  in  the  brackets  {}  does  not  involve  the  0's  or  the  0's, 
we  can  perform  the  integrals  over  these  variables,  which  result  in 
Fourier  transforms  of  the  original  intensity  distributions.  Let  the 
Fourier  transform  of  0(0,0)  be  0  (kz,kw).  Then  we  can  rewrite  Eq.  (A-5): 

a  2  ( £+m-l )  21  2m  I  ~  12  iki  .  ..} 

=  (2n) J  Xdu1  dz  du^  dv^  dw  dv^  (k  f)  z  w  |0(kz,kw)|  e  , 

(A-6) 

where  the  term  in  brackets  (...)  is  the  same  as  in  the  preceding  equation. 
Except  for  the  exponential,  the  integrand  in  Eq.  (A-6)  is  positive- 
definite.  Therefore  the  integral  will  reach  its  maximum  value  when  the 
exponential  is  identically  equal  to  one — i.e.,  the  term  in  brackets 
{...}  =0.  This  term  is  identically  zero  only  if  6(u,v)  is  zero,  or 
at  most  a  linear  function  of  u  and  v: 


fi(u,v)  =  a  +  bu  +  cv 


(A-7) 


where  a,  b,  and  c  are  constants.  A  linear  variation  of  6(u,v)  results 

in  a  shift  in  the  image  position,  but  in  no  image  distortion.  Therefore 

we  have  proven  that  maximizing  S  restores  all  but  the  original  position 

4 

of  the  undiscorted  image. 
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2.  Proof  that  Sg  Reaches  an  Absolute  Maximum  for  an  Unresolved  Star 

In  Table  1  of  the  main  text  we  defined  S  : 

5 

S  =  /  dx  dy  I  (x,y)  n  ^  2  ,  integer  .  (A-8) 

o 


For  an  unresolved  star,  0(0  ,<£)  is  a  Dirac  delta  function,  which  for 
convenience  we  put  at  0  =  0  and  <p  =  0.  Then,  going  through  a  derivation 
similar  to  that  given  for  S^,  we  can  rewrite  Eq.  (A-8): 


S 

5 


/2jrf\  n  r  n  n  n-1  n-1  ikj  .  .  .  } 

I — ^ —  )  0(0,0)/dudvd  z  d  we 


( A-9) 


where  the  brackets  {...]  = 


n-1 


Z[  6  ( u  . ,  v  . ) 

l  1 

j=l  J  J 


6 (u ,-z  , v  -w  )  |  +  6 (u  , v  ) 
J  J  J  J  n  n 


and  contain  the  effect  of  the  atmospheric  perturbation.  As  in  the 
preceding  proof,  S  will  have  its  absolute  maximum  only  when  {...}  =  0, 
which  again  implies  that  6(u,v)  is  at  most  a  linear  function  of  its 
arguments — i.e.,  when  the  atmospheric  distortion  has  been  removed. 

3,  Maximum  Theorem  lor  S3  and  Sg  (Image-Defect  Functions) 


S  and  S  are  examples  of  sharpness  definitions  in  which  the 
*3  o 

observed  intensity  I(x,y)  is  compared  with  a  standard  of  reference. 
If  we  have  advance  knowledge  of  the  undistorted  intensity  distributi 
IQ(x,y) .  then  we  can  attempt  to  maximize  the  image-defect  function  S 


on 

8  ’ 


S  =  -  /dx  dy  I  (x,y)  -  I  (x,y) 
o  U 


( A- 1 0 ) 
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This  function  obviously  reaches  its  maximum  value  (zero)  only  when 
I ( x , y )  exactly  matches  IQ(x,y).  If  we  expand  the  square  we  get 

9  -  2 

S  =  -  /  dx  dy  l“(x,y)  +  2  /dx  dy  I^x.y)  l(x,y)  -  /  dx  dy  l^x.y) 

(A-ll) 

The  first  term  in  this  equation  is  -S  ;  the  second  term  is  2S^  ior  the 

particular  case  when  M(x,y)  =  I^(x,y);  and  the  third  term  is  a  constant. 

Since  S  is  maximized  at  the  restored  image  and  the  S  is  also  maximized, 

8  1 

it  follows  that  must  also  reach  its  absolute  maximum  if  Eq.  (A-10)  is 
to  be  correct . 

We  have  not  made  analytic  studies  of  what  happens  if  the  function 

M(x,y)  in  does  not  match  l^x.y),  except  in  the  case  where  M  is  a 

Uirac  delta  function,  in  which  case  reduces  to  S^,  and  in  the  case 

where  M(x,y)  =  -x2  -  y2  (for  which  S3  is  equal  to  S6  of  Table  1),  where 

24 

a  maximum  proof  has  been  obtained.  Our  computer  simulations  indicate 
that  M(x,y)  can  depart  from  I^(x,y)  in  substantial  ways  before  our  iteia- 
tion  technique  fails  to  converge  on  a  true  image.  All  the  defect  functions 
have  the  feature  that  they  automatically  center  the  image  as  well  as 
sharpen  it. 
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